Ratchet transport of overdamped particles in superimposed driven lattices
Huang Shu-Na, Zhu Wei-Jing, Huang Xiao-Qun, Ai Bao-Quan, Li Feng-Guo
Guangdong Provincial Key Laboratory of Quantum Engineering and Quantum Materials, School of Physics and Telecommunication Engineering, South China Normal University, Guangzhou 510006, China

 

† Corresponding author. E-mail: aibq@scnu.edu.cn

Abstract
Abstract

Ratchet transport of overdamped particles is investigated in superimposed driven lattices using Langevin dynamics simulations. It is found that noise can strongly affect the transport of the particles. When lattices driving dominates the transport, the noise acts as a disturbance of the directed transport and slows down the average velocity of the particles. When the driving phase has less impact on particle transport, Gaussian white noise can play a positive role. By simply modulating these two parameters, we can control efficiency and the direction of the directed currents.

1. Introduction

In recent years, the rachet effect in driven lattice potentials—which can evoke directed particle transport from an unbiased nonequilibrium system—has been a popular research topic.[18] The rachet effect is the phenomenon of driving thermal Brownian particles into directed currents by breaking certain space–time symmetries in the absence of a net force. Broadly speaking, there are four types of rachet models: rocking ratchets,[9,10] flashing ratchets,[1113] correlation ratchets,[1417] and entropic ratchets.[1821]

Particles can be pumped in periodic asymmetric potential by switching the sawtooth dielectric potential on and off.[22] Further study finds current reversal occurring in deterministic ratchets, which is related to bifurcation from a chaotic and periodic regime with the change of an exponent and a control parameter.[23] A symmetric periodic potential and a zero-mean friction force produced by two counterpropagating lattice beams can also break the temporal symmetry of the system and induce directed motion, which has been demonstrated experimentally.[24] It has been theoretically and experimentally demonstrated that the interaction between deterministic driving and fluctuations produces the resonance which is related to driving frequency. The direction of the transport is determined by the driving frequency.[25] Based on the mechanism of ratchet effect, optical lattices can be used to trap colloidal Brownian particles and induce systematic transport in experiments.[2628] Brownian particles can be transported in a two-dimensional (2D) coupled flashing ratchet driven by Gaussian noise.[12] In the anisotropic structures, white thermal noise together with periodic driving force can generate a non-zero macroscopic velocity.[9] Rectification phenomenon specific to high-dimensional rocking ratchets is determined theoretically and experimentally.[29,30] As for deterministic optical rocking ratchets, the direction of current can be changed by the asymmetry of the potential.[31] Experiments show that dissipation can induce currents for cold atoms in a space–time symmetric potential.[32] In a 2D ac-driven lattice, changing the structure of the lattice in the direction perpendicular to the applied driving force can induce current reversals.[33] Another study finds out using a periodically shaken 2D dissipative lattice can separate particles in specific directions by their physical properties.[34]

Recently, Mukhopadhyay and coworkers superimposed a symmetric oscillating lattice as “substrate lattice” on a second oscillating lattice as “carrier lattice” to manipulate particle transport in real time. By switching the carrier lattice on (off), directed current can be accelerated (freezed). Directed current can be slowed down by changing the value of a phase difference.[35] However, noise is not considered in this work. As we know, noise plays an important role in rachet transport. In this paper, we study the directed transport of overdamped particles in the superimposed driven lattices and focus on finding how noise affects the rectification.

2. Model and methods

We consider N noninteracting overdamped particles moving in a periodic potential as shown in Fig. 1. Here, Vs represents substrate lattice and Vc represents carrier lattice.[35] The Vs and Vc are

where k is the wave number. The oscillation amplitude and frequency of carrier lattice are twice as large as those of substrate lattice. Each particle obeys the overdamped Langevin equation
where x is the position of the particle and γ is the friction coefficient; D0 is the noise intensity; and the Gaussian white noise is ξ(t) with zero mean and unit variance. The average velocity in the asymptotic long-time regime can be obtained from the formula,

Fig. 1. Schematic illustration of particles moving in superimposed driven lattices called substrate lattice and carrier lattice.
3. Results and discussion

In the simulation, the integration time step is set as and the total integration time is set as 104. Unless otherwise noted, our simulations are performed under the conditions of the parameter sets: , k=1.0, , and γ=1.0. Our simulation results do not depend on time step and integration time when and . We vary D0, d, ω, and ϕ and calculate the average velocity V of passive particles. The simulation results are shown in the following Figs. 27.

Fig. 2. Average velocity V as a function of noise intensity D0 for different ϕ at d=π and ω=1.0.

Figure 2 shows the average velocity V as a function of translational diffusion coefficient D0 for different ϕ. When ϕ=π, the directed motion cannot be observed. For ϕ=π/2 and 3π/2, the directed current is most significant and the average velocity V decreases monotonically when D0 increases. In these cases, lattices driving dominates the transport, and transports the particles efficiently. With the increase of D0, particles diffuse faster and faster randomly, which means the noise acts as a disturbance and is against the directed transport. For ϕ=π/2 and 3π/2, the average forces provided by lattices are completely opposite in a time period 2π, which indicates that the time cumulative effects of lattices are completely opposite and induces opposite currents. For ϕ=π/4 and 5π/4, there exists an optimal value of D0 ( ) at which V reaches its maximal. For ϕ=π/4 and 5π/4, the rachet effect becomes faint, and the particles diffuse faster and are easier to escape from the well with the D0 gradually increasing from 0. When particles can just enough pass through the well, V reaches to a peak which is concerned with the oscillation characteristics of the lattices. When , the rachet effect disappears and the particle starts the chaotic motion.

Figure 3 shows the average velocity V as a function of d with or without white noise. As shown in Fig. 3, when , oscillation amplitude is so small that lattices driving disappears. When , the oscillation amplitude is too large, thus there is not enough time for particles to diffuse away. From Fig. 2, we can find out when lattices driving dominates the transport, noise acts as a disturbance of the directed transport and slows down the average velocity of the particles, and shown in Fig. 3. For example, when ϕ=π/2 and d = 3.7, V takes its maximum up to 0.5 without noise (see Fig. 3(a)), but it only takes 0.4 for overdamped particle (see Fig. 3(b)). When the rachet effect becomes faint, particles diffuse faster with noise. Therefore, noise facilitates the directed transport such as the curve at ϕ=π/4 and d=4.7. In addition, when ϕ=π, the system is time-reversal symmetry with or without noise. Only the Gaussian white noise cannot induce a particle current. In Fig. 3(a), when d is appropriate, the oscillating lattices show the feature of resonance and cause the jump of data. The resonance condition is determined by d, ω, and ϕ, as shown in Figs. 3,5, and 6. In particular, we further verify that particle number N has almost no impact on our results as shown in Fig. 4.

Fig. 3. Average velocity V as a function of oscillation amplitude d for different ϕ. (a) Without white noise. (b) With white noise. The other parameters are D0 = 0.1 and ω = 1.0.
Fig. 4. Average velocity V as a function of particle number N at ϕ=π/2, d=π, D0 = 0.1, and ω=1.0.
Fig. 5. Average velocity V as a function of oscillation frequency ω for different ϕ. (a) Without white noise. (b) With white noise. The other parameters are D0 = 0.1 and d = π.
Fig. 6. Average velocity V as a function of driving phase ϕ for different d. (a) Without white noise. (b) With white noise. The other parameters are D0 = 0.1 and ω = 1.0.
Fig. 7. (a) Contour plots of the average velocity V as a function of d and ϕ with white noise at D0 = 0.1 and ω = 1.0. (b) Contour plots of the average velocity V as a function of ω and ϕ with white noise at D0 = 0.1 and d = π.

Figure 5 displays the average velocity V as a function of ω with or without white noise. When , the lattices stop oscillating and the particle pump is not allowed. When , the lattices oscillate so fast that the particles cannot diffuse away. In Fig. 5(a), particles are transported to positive direction only by lattices driving at appropriate ω, while in Fig. 5(b), there also exists noise driving, thus these two kinds of driving will compete. For ϕ=3π/2, particles take the maximal absolute average velocity when oscillation frequency ω is nearly 1 with or without noise, which means lattices driving dominates the transport and noise disturbs the directed motion. Especially, when ω is much smaller (or larger) than 1, the effect of lattices driving recedes and noise driving can facilitate the transport. For ϕ=π/4, we can arrive at the same conclusions. Therefore, the value range of ω allowing for the directed current is broadened with the white noise (see Fig. 5(b)), and ω can be the switch of the current reversal.

Figure 6 describes the average velocity V as a function of ϕ with or without white noise. As the figures show, all curves are antisymmetric about π. Therefore, ϕ can be the switch to control the direction of the current. When d=3π/4, π, and 5π/4, oscillation amplitude reaches a appropriate value and lattices driving dominates the transport. In these cases, noise seems to slow down the directed currents. However, when d=π/2 and 3π/2, no current occurs without the noise. The oscillation amplitude is out of limits, thus white noise plays an important role on directed transport when lattices driving is ineffective.

To study more details of the dependence on these parameters, we calculate the average velocity V as a function of d and ϕ in Fig. 7(a). When d=d1, average velocity V is non-negative for and non-positive for , thus ϕ can switch the direction of the current. For , it changes its direction twice when d is increased from zero. Figure 7(b) shows the contour plots of average velocity V as a function of ω and ϕ. For ω=ω1, the average velocity is positive when and negative when . For and , it changes its direction once when ω is increased from zero.

4. Conclusion

To summarize, we study the transport of noninteracting particles moving in superimposed driven lattices. We calculate the dependence of noise intensity D0, oscillation amplitude d, oscillation frequency ω and driving phase ϕ on the directed currents. Base on Mukhopadhyay and coworkers’ study, we further explore three aspects in this work.

(i) We consider noise in this model and focus on the difference between transporting particles in superimposed driven lattices without and with noise.

(ii) We find that noise plays an important role in superimposed driven lattices and get some interesting result.

(iii) We also provide contour plots of the average velocity and show the dependence of these parameters collectively and directly.

It is found that the influence of noise driving on the rectification is relative to the driving lattices. When lattices driving dominates the transport, the noise acts as disturbance of the directed transport and slows down the average velocity of the particles. When the driving phase has less impact on particle transport, Gaussian white noise can play a positive role. By simply modulating these two parameters, we can control efficiency and the direction of the directed currents.

Reference
1 Salger T Kling S Hecking T Geckeler C Morales-Molina L Weit M 2009 Science 326 1241
2 Denisov S Flach S Hänggi P 2006 Europhys. Lett. 74 588
3 Brown M Renzoni F 2008 Phys. Rev. 77 033405
4 Gommers R Lebedev V Brown M Renzoni F 2008 Phys. Rev. Lett. 100 040603
5 Reimann P 2002 Phys. Rep. 361 57
6 Hänggi P Marchesoni F 2009 Rev. Mod. Phys. 81 387
7 Dittrich T Dubeibe F L 2015 Phys. Rev. Lett. 114 094101
8 Cubero D Renzoni F 2016 Phys. Rev. Lett. 116 010602
9 Bartussek R H?nggi P Kissner J G 1994 Europhys. Lett. 28 459
10 Magnasco M O 1993 Phys. Rev. Lett. 71 1477
11 Reimann P 1997 Phys. Rep. 290 149
12 Bao J D Zhuo Y Z 1998 Phys. Lett. 239 228
13 Reimann P Bartussek R Haussler R Hänggi P 1996 Phys. Lett. 215 26
14 Doering C R Horsthemke W Riordan J 1994 Phys. Rev. Lett. 72 2984
15 Bartussek R Reimann P Hänggi P 1996 Phys. Rev. Lett. 76 1166
16 Büttiker M 1987 Zeitschrift für Physik B-Condens. Matter 68 161
17 Usmani O Lutz E Büttiker M 2002 Phys. Rev. 66 021111
18 Ai B Q Liu L G 2006 Phys. Rev. 74 051114
19 Ai B Q 2009 Phys. Rev. 80 011113
20 Marchesoni F Savel’ev S 2009 Phys. Rev. 80 011120
21 Ai B Q 2009 J. Chem. Phys. 131 054111
22 Rousselet J Salome L Ajdari A Prostt J 1994 Nature 370 446
23 Mateos J L 2000 Phys. Rev. Lett. 84 258
24 Schiavoni M Sanchez-Palencia L Renzoni F Grynberg G 2003 Phys. Rev. Lett. 90 094101
25 Gommers R Douglas P Bergamini S Goonasekera M Jones P H Renzoni F 2005 Phys. Rev. Lett. 94 143001
26 Faucheux L P Bourdieu L S Kaplan P D Libchaber A J 1995 Phys. Rev. Lett. 74 1504
27 Lee S H Ladavac K Polin M Grier D G 2005 Phys. Rev. Lett. 94 110601
28 Lopez B J Kuwada N J Craig E M Long B R Linke H 2008 Phys. Rev. Lett. 101 220601
29 Denisov S Zolotaryuk Y Flach S Yevtushenko O 2008 Phys. Rev. Lett. 100 224102
30 Lebedev V Renzoni F 2009 Phys. Rev. 80 023422
31 Arzola A V Volke-Sep’ulveda K Mateos J L 2011 Phys. Rev. Lett. 106 168104
32 Gommers R Bergamini S Renzoni F 2005 Phys. Rev. Lett. 95 073003
33 Mukhopadhyay A K Xie T Liebchen B Schmelcher P 2018 Phys. Rev. 97 050202
34 Mukhopadhyay A K Liebchen B Schmelcher P 2018 Phys. Rev. Lett. 120 218002
35 Mukhopadhyay A K Liebchen B Wulf T Schmelcher P 2016 Phys. Rev. 93 052219